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UNIT I 

PROGRAM ABSTRACT AND HISTORY OF USE 


A, ABSTRACT 

Program THNSKN computes the mean heating rate at a 
number of thermocouple locations on the surface of 
thin skin transient heating rate models. Experi- 
mental data, calibrated in engineering units, is 
obtained from ABEL tape output of the EHTS03 data 
reduction computer program. The program output con- 
sists of time history tabulations in an array by 
thermocouple locations for each time sample. 

B. PROGRAM HISTORY 

The original THNSKN program was written in 1969. 

The current version was written in December 1974 as 
requested by SCR83-15-30. 




UNIT II 

PROGRAM USERS INFORMATION 


A. IDENTIFICATION 
Title : 

Author : 

Date : 

Installation : 
Authorization : 


Surface Heating Rate of Thin 
Skin Models (THNSKN) 

Jack D. McBryde 

December 1974 

JSC, Houston, Texas 

Contract NAS 9-12200 


Source Language: FORTRAN V 

Computer Configuration: UNIVAC 1108 or 1110, EXEC 8 

B. DESCRIPTION OF PROGRAM REQUEST 
See attachment 8. 

C. USAGE 

1. Input Requirements and Data Descriptions: 

The THNSKN input tape is produced by the program 
EHTS03. It is a FORTRAN-written binary tape 
containing time in floating-point seconds fol- 
lowed by data in floating point for N channels. 
The number of channels CN) is fixed for a given 
tape (N £ 259) . 

See attachment 6 for sample input. 

2. Lead Card Setup: 

See attachment 1. 
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3 . Deck Setup : 


See attachment 2. 

4. Tape Assignments: 

See attachment 3. 

5. Restrictions: 

The following list is a summary of the restrictions 

• The input tape must have the format specified 
in section C.l. 

• The array of thermocouples must have no more 
than 50 rows or columns, and one row or column 
must contain at least three thermocouples. 

• Each specified time slice must contain at least 
three input records. 

• The maximum number of points for thermocouples, 
which may be used for least squares heat loss 
computations, is 50. 

6. Diagnostic, Normal-Console, and Standard-Output 

Messages: 

See attachment 9. 

7. Labels, Save-Tapes, and Output Disposition: 

There is printer output only. 

8. Output Formats: 

See attachment 7 . 


3 


9. Nonstandard System Requirements: 

There are none. 

10. Running Time: 

THNSKN processes nine thermocouples for 40 time 
points in less than 1/2 minute. 

D. ANALYSIS 

1. Description of the Relationship of Variables to 
the Test Item: 

The test item may be considered to be a rectan-: 
gular plate having, at most, two thicknesses as 
shown : 


y ► 



• The origin of the x- , y-coordinate system is the 
upper left corner. 
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• The variables and d 2 are the two thick- 

nesses. If the test item has only one thickness, 
it is represented by d . 

• The variable T. . is the temperature at a 

3 

thermocouple in degrees Fahrenheit. 

• The variables T^ k and T i k+1 are tempera- 
tures adjacent to the thickness difference 
junction. 

• The variable Y. , is the y-coordinate of the 

1 » d 

junction at row i in inches. 

• Each row or column is colinear in the X-, 

Y-plane so that the coordinates of all the 
thermocouples (X, Y) can be expressed using 
i-values of X and j -values of Y . Both X 
and Y are in inches . 

2. Equations: 

2.1 Mean temperatures: 

The mean temperature T for each thermocouple is 
computed by 


N 

^ = N T £ (1) 

£=1 

where N is the number of temperatures' read in 
the time interval processed. 
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2.2 Mean heat flux without conduction correction: 


The mean heat flux without conduction correction 
q is computed by 


L i',3 


pCp( 

24(t^ 


N - 1 / t N 
t x ) N - 1 \ i > j 


t N-1 

i,3 


•t} 


i»3 


- T 


2 •) 
1 » 3 / 


( 2 ) 


where 

t^ 


'N 


. , T? T 1 ?'*, T 1 ? . 
1,3 1,3 i,3 i,3 


P 

C, 


= the time for the first 

scan of data in the interval 
in seconds. 

= . the time for the last scan 
of data in the interval in 
seconds . 

= the first, second, next to 
the last, and last tempera- 
tures of the thermocouple in 
row i and column j . 

= mass density in lb/ft . 

= specific heat in BTU/lb°F. 

= thickness of test material 
in inches . 


If the test item has two thicknesses, then d = d^ , 
for j <. k (see the figure) and d = d^ for 
j > k . 
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2.3 Mean lateral heat loss in the X-direction: 

The mean lateral heat loss in the X-direction q 
is computed by 


1 » 3 


24Kd 

X i+1 " X i- 1 


T . 


- T 


T. 


X. - X. , 
l l-l 


i ~ 1 >j . 1+1 >3 


X i + 1 


T. 

1 > 3 
X . 


C3) 


where 

K = thermal conductivity in BTU/ft sec°F. 
d = d^ for j <. k (see the figure) . 

d = f° r j > k • 

2.4 Mean lateral heat loss in the Y-direction: 

The mean lateral heat loss in the Y-direction 
may be computed by one of three equations, 
depending on whether the test material has one or 
two thicknesses and the position of the thermo- 
couple relative to the junction between the 
thicknesses . 


Case 1 : The variable q^ when the material has 

one thickness or where it has two and all three 
thermocouples are the same side of the junction. 


q i. j 


24Kd 


T. 


3+1 




JL»JL 

Y j 


- T 


i z 1 


Y. , 
3-1 


hllL 

y 3 + 1 


- T. 




Y. 

3 


(4) 
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where 

d = dj when j + 1 <. k , 
d = d 2 when j - 1 > k . 

Case 2 : The variable when the material has 

two thicknesses and j ■ = k (see the figure) . 



24Kd 1 

Y k.i - Vi 



- T 


i ,k- 1 


k-1 


T r 


2^ T i,k+l 




‘ T i k 3 
^k + 'l'-"*d> d 


3 


( 5 ) 


where 

= the Y-coordinate of the function at row i 
(see the figure). 

Case 3 : The variable q^ when the material has 

two thicknesses and j = k + 1 . 


q 


y 

i J 


24Kd 2 
Y k + 2 - Y k 


d l< T i,k + l - T i,k^ 

t Y d - V k )ff 2 ♦ - -V^d^ 


T 

i,k+2 

— ¥ 7 — 

*k+2 


T i ,k+l 


( 6 ) 


2.5 


3. 


4. 


5 . 


Mean heating rate: 

The mean heating rate q is computed from 
equations 2, 3, and 4, 5, or 6. 






(7) 


Special Treatment When the Array Contains Bad 
Thermocouples: O 

The variables q, q, q X and q^ are not computed 
for bad thermocouples. These values are set to 
zero . 

In computing q x for a column containing one or 
more bad thermocouples, temperatures and 
x-coordinates of three successive good (but not 
necessarily adjacent) thermocouples are used. 

The same procedure is used in computing q 7 . 

The Variables q X and q^ at Row or Column End 
Points: 

X v 

The conduction correction (q or oj ) used at the 
end points (first or last good thermocouple) are 
the same as those computed for the adjacent 
thermocouple in the row or column. 

Linear Least Squares Heating Rate and Percent 
Standard Error : 

' • 

The heating rate at each thermocouple Q is com- 
puted from the slope M by 


9 


Q - 


pCpd 

12 


M 


( 8 ) 


The slope M is computed by 


M 


n£Vi - St. Et, 

n I>i - (EV 2 


C9) 


where 

n = the number of scans, 
t = time. 

T = temperature. 


The percent standard error E is computed by 


s n (M) 

E = -Sjj— x 100 


(10) 


S n (M) = na n CT)J(n - 2)[n£t 2 - (^tp 2 (11) 


a (T) = - 
n v J n 


"Et 2 - tEip 


2 . ^EVi - E^Ey 

n E*- - c£y 2 


2 .!/ 2 


( 12 ) 
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If the time slice specified contains more than 
n scans, the last time and corresponding 
temperatures are used as the first time and 
temperatures' for the next n scans. 

E. FLOW CHART 

See attachment 4. 

F. STORAGE REQUIREMENTS 

Forty-three thousand decimal locations are required. 

G. LIST OF SUBROUTINES 

None. ■ 

H. LOCATION OF SOURCE AND BINARY DECKS 

The Data Processing Systems Department of Lockheed 
Electronics Company, Inc., maintains all source and 
binary decks for this program. 

I. PROGRAM LISTING 
See attachment 5. 

J. CORRESPONDENCE 
See attachment 8. 
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ATTACHMENT 1 
LEAD CARD SETUP 






■ ;i > 

LEAD CARD SETUP 


CARD NO._I 

JOB PAGE 1 OFJJ. 

NAME T ^ S ,^ PROGRAMMER McBr yd e DATE 12/ 7 4 


FIELD 

ID 

CARD 

COLUMNS 

FORMAT 

SYMBOLIC 

NAME 

IDENTIFICATION 

1 

1-72 

12A6 

TITLE1 

First title card. 




{ 

LEAD CARD SETUP 

• 

CARn NO. 2 



JOB 


PAGE 2 ft F 10 

name thnskn 

PROGRAMMER McBryde 

DATE .127-74 


FIELD 

ID 

CARD 

COLUMNS 

FORMAT 

SYMBOLIC 

NAME 

IDENTIFICATION 

1 

1-72 

12A6 

TITLE2 

Second title card. 


























































































- " . * ■ ' 
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CONSTANTS 


LEAD CARD SETUP 


CARD NO._i 

JOB. PAGE 3 of ...10 

name THNSKN PROGRAMMER McBryde . DATE 1 j/_L 4 .. 


field! 

id! 

CARD 

COLUMNS 

FORMAT 

SYMBOLIC 

NAME 

IDENTIFICATION 

H 

1-10 

F10.3 

RHO 

Mass density (p) of test material in 

II 




lbm/ft^. 

m 

11-20 

F10 . 3 

CP 

Specific heat (Cp) of test material in 

■ 




BTU/lbm °F . 

m 

21-30 

mi 

BIGK 

Thermal conductivity ex’) of test 

■ 




material in BTU/ft sec °F. 

H 

31-40 

F10.3 

DD(1) 

Test material thickness (d^) in inches. 

| 




If there is only one thickness in the 

■ 




sample, the next field will be left 

■ 




blank. 

5 

41-50 

F10.3 

DD (2) 

The second test material thickness (D ? ) . 

■ 




If a second thickness is supplied, 

■ 




card 5 must appear in the deck. 

6 

51-53 

13 

NX 

Number* of x-coordinates : 

■ 




1 < NV < 50 . 

7 

54-56 

13 

NY 

Number* of y-coordinates: 

■ 



. 

■ 

1 < NY < 50 . 

■ 





■ 





■ 

Pfg| 




■ 





■ 

■1 





Comment: *A minimum of a 1x3 or 3x1 array of x, y-coordinates is 

required. 
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MISCELLANEOUS RUN PARAMETERS 


LEAD CARD SETUP 

CARD KO._* 

JOB 

NAME THNSKN PROGRAMMER McBl 7 de 


FIELD 

ID 

CARD 

COLUMNS 

FORMAT 

SYMBOLIC 

NAME 

IDENTIFICATION 

1 

1-5 

15 

NCHS 

Number of data channels on input tape*. 

2 

6-10 

15 

NCC 

Command channel number*. 

3 

11-15 

15 

IOPT1 

• ■ ■ 

0 — Use card 9 to specify start-stop 





times . 





1 — Use card 9A to specify start-stop 





command event levels. 

4 

16-20 

15 

NPAR 

The maximum number of parity error 





records per time slice which will be 





permitted. If blank, 2 will be used. 





Parity error records are skipped. 

5 

21-25 

15 

NOTH 

If NOTH is blank or zero time, tempera- 





ture data for each scan in the time 





slice will be printed. 











.. 




































Comment: *The input tape record has a time word preceding the 

data channels. These numbers do not include the time 
word. 

1-4 


PAGE_i OF_i2 

DATE 12/-Z4 



















Y- COORDINATES OF THICKNESS DIFFERENCE JUNCTION 


CARD NO._L 

JOB 

NAME 


CARD BftBMAT SYMBOLIC 
COLUMNS PORMA1 NAME 


LEAD CARD SETUP 



PAGE § OF ,.1-Q 

PROGRAMMER McBryde _.dkje1I/2±- 


IDENTIFICATION 


NOTE: This card will appear in the 


deck only if DD(2) on card 3 


is used. 




11-20 

F10.3 

3 

21-30 

F10 . 3 

m 

31-40 

F10.3 

5 

41-50 

F10.3 


^-coordinate of the junction 


between test material thicknesses when 


the x-coordinate is xxfl' 



Comment: Use as many cards as required to define NX values 

1-5 


























X- COORDINATES 


LEAD CARD SETUP 


FIELD CARD format SYMBOLIC 
ID COLUMNS " A NA mE 


AME THNSKN programmer McBryde 


IDENTIFICATION 


1-10 F10.3 XX (1) Ix-axis location of first thermocouple 


row in inches . 


2 11-20 F10.3 XX(2) X-axis location o£ second thermocouple 


row. 


X-axis location of third thermocouple 


3 

21-30 

F10.3 

XX 


31-40 F10.3 XX (4) 


41-50 F10.3 XX (5) 


X-axis location of fourth thermocouple 


row. 


X-axis location of fifth thermocouple 


row. 



Comment: Enter as many cards as needed to define NX locations 

(NX < 50). 



































Y- COORDINATES 


CARD NO.. 
JOB 

7 

LEAD CARD SETUP 

PAGE 7 

o 

i—i 

u. 

C 

NAME 

THNSKN 

PROGRAMMER McBryde 

DATE 

12/74 


FIELD 

ID 

CARD 

COLUMNS 

FORMAT 

SYMBOLIC 

NAME 

IDENTIFICATION 

1 

1-10 

F10.3 

YY (1) 

Y-axis location of first . thermocouple 





column in inches. 

2 

11-20 

F10.3 

YY (2) 

Y-axis location of second thermocouple 





column. 

3 

21-30 

F10.3 

YY (3) 

Y-axis location of third thermocouple 





column. 

4 

31-40 

F10.3 

YY (4) 

Y-axis location of fourth thermocouple 





column. 

5 

41-50 

F10.3 

YY (51 

Y-axis location of fifth thermocouple 





column. 






























































Comment: Enter as many cards as required to define NY locations 

(NY < 50). 
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THERMOCOUPLE CHANNEL NUMBERS 


LEAD CARD SETUP 


PAGE_? OF 10 


FIELD CARD 
ID COLUMNS 


THNSKN 


FORMAT 


13 


SYMBOLIC 

NAME 



PROGRAMMER McBryde 


IDENTIFICATION 


Row number (X-direction) o£ 


thermocouple. 


DATE 1 .L // 7 - 4 .- 


99 

13 

IY 

Column number (Y-direction) 



thermocouple. 


Channel number for temperature of 


thermocouple located at IX, IY, 


y Second thermocouple. 



Third thermocouple. 



Fourth thermocouple 



Fifth thermocouple. 



Comment : 


Enter as many cards as needed to define NX * NY thermo- 
couple locations. Bad thermocouples should be indicated 
with NTCH (IX, IY) = 0 . 




















































START -STOP TIMES 


CARD NO.. 
JOB 

9 

LEAD CARD SETUP 

PAGE 9 OF 10 

NAME 

THNSKN 

programmer McBryde 

DATE 12/74 


FIELD 

ID 

CARD 

COLUMNS 

FORMAT 

SYMBOLIC 

NAME 

IDENTIFICATION 

1 

1-10 

F 1 0 . 3 

TSTRT 

Computation start time. 

2 

11-20 

F10.3 

TSTOP 

Computation stop time. 

3 

21-30 

F10.3 

STATME 

Start time for least squares curve fit. 





No data will be processed outside of 





the interval from TSTRT to TSTOP. 

4 

31-35 

15 

NTME 

The number of points to be curve fit. 





If NTME < 2 , no curve fitting will 





be done. 






■ 




NOTE: This card will be used if IOPT1 

■ 




on card 4 is zero or blank. Card 9A 

■ 




HHHHHHMMHWi 

pi 




in the deck. As many start- stop time 

■ 




cards may be used as desired but the 

■ 

Hi 




■ 




start time for the current time slice 





must be >. the stop time for the 

■ 




previous time slice. 

■ 





■ 

(Hi 




■ 

hi 




■ 

Hi 

















COMMAND EVENT START- STOP 


CARD HO._JA 

JOB 

NAME THNSKN 


LEAD CARD SETUP 


PROGRAMMER McBryde 


PAGE 10 OF 10 
DATE 12 / 74 - 



CARD 

COLUMNS 

FORMAT 

SYMBOLIC 

NAME 

IDENTIFICATION 

m 

1-10 

F10.0 

CESTRT 

Command event start level. 

2 

11-20 

F10.0 

CESTOP 

Command event stop level. 

3 

21-30 

F10.3 

mmmm 

Start time for least squares curve fit. 


31-35 15 


NTME 



No curve fit computations will be done 


prior to CESTRT or after CESTOP. 


The number of points to be curve fit. 


If NTME <_ 2 , no curve fitting will 


be done. 


NOTE: Card 9A will appear in the deck 


only if IOPT1 on card 4 is 1. Card 9 


ear in the decks. As man 


event start-sto 


required, but intervals requested must 
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ATTACHMENT 2 
DECK SETUP 



DECK SETUP FOR PROfiWAM -THNSKN (EXEC VIII 


(Bock of dock) 



One or more card 9A’s if I0PT1 
on card 4 is 1, 


One or more card 9*s if IOPT1 
on card 4 is 0. 


CARDS 6-8 


CARD 5 [Only if DD(2) on card 3 > 0] 


THNSKN CARDS 1 



8XQT . THNSKN/ ABS 


8MAP.S . THNSKN , THNSKN/ ABS 


/ 7 

1 8PREP 

7 

8C0PIN PGMTP ,TPF$ 


8ASG,T 8.8C.XXXXXX (Input tape number) 


8ASG,T PGMTP, 8C,V02679 


8MSG.N TWO TAPE. DRIVES 


7 

8 RUN 



(Front of dock) 


( ) 


ATTACHMENT 3 


RUN REQUEST FORM 


INSTRUCTIONS FOR CENTRAL COMPUTER COMPLEX COMPUTER RUNS 

mm (DO MOT FILL IN SMAOCO AAC AS) 

MOOMMIIII COMMtN TSl 

EXEC VIII 



J!.C tl JHM SfXA Rl t MAT /'.*» 


a u s. commnur wiuting oirici u/4— 7/1 007 u n 


3-1 



















ATTACHMENT 4 
FLOW CHART 



THNSKN 

















ATTACHMENT 5 
PROGRAM LISTING 


| 


a F_0« 1 t_s »THN5KN i •THNSk'N 

FOR SE I X-l 2/20/71-03: 021 19 (,0> 

PROGRAM " ~ " ’ 

STORAGE^ used: CODE(l) 003717: DATA(O) 073316: blank C0MM0N<2> 000000 



0017 NERR9S 

0020 NERR3S 


STORAGE ASSIGNMENT 1 BLOCK » TYPE, RELATIVE LOCATION, NAME) 


oooj oooiso 10 L Ero^oo _ 0723 3 s iooof 

Cn '"0000 ~ 0 7 2 3 ** 3 iooif" ~ noon" r r 7 T 3 '*fS — ioOSF 

1 0001 001711 I P22G 0001 001750 10276 

H* • 0001 00047? 1 10 L* oo'ol 'CTCJTTdY THTTg" 

0001 002151 11316 0001 002155 Il37« 

0001 002107 1 2 1 1 G 0001 SoTTl 4 7 251 G 

0001 003017 1 3 7 7 G 0001 000732 1H0L 

oooi~~ bo3"f57 " 'i**53g ' odbl 003171 - n*TG‘ 

0001 00323* 1S06G 0001 

nooi "" oo33Trr*isi2G ooor 

0001 003515 I A 7 7 G 000) 


00 00 

0000 

0001 


07 2336 1 00 1 F 
072350 1 007F 

001771 1011S 

— rmry 


0007 O02TT7 
002 t 75 
002116 


0001 

oOoi 


77770" 

1 1536 


0001 

003632 

20 OL 

oooo 

072561 

2001F 

oooo 

' 0? 7663" 2no?F ’ 

oooo 

072736 

20 1 IF 

o 0 oft” 

'"O7503T 3o19f 

0001 

001056 

2?0t. 

0001 

oorn^ 2 iol 

0001 

001 126 

2 AOL 

0001 

00 1 176 

30 OL 

0001 

001363 

365 L 

' 5000 

073735 

1OT2f“ 

0001 

00061 1 

107G 

• "0007 507070’ 

172T5 

0001 

0023SS 

5 1 OL 

0001 

002173 

S50C " 

OOP l 

0025.57 

575L 

0001 

(102717 

600 L ' 

nooi 

003131 

610L 


0000 

°?°5L_ 

oooo 

0000 
0000 

0001 

' ctoot~' 
000 1 
000 1 ' 
OOOI 

0000 

0001 

-ootn — 

oun v 
bo'oi 
0001 
0001 ' 
0001 


000732 

OoJTTr 

003215 
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00919 

19 9* 

00915 

200* 

00916 

201* 

00917 

202* 

00920 

203* 

00921 

209* 


IF( I0PT1 .NE.O) IOPT1 « 1 
' W R I f E • 6 » 2 0 0 2 ) N C H S • N C C , I OPT 1 ,NPAR .NOTH 
IF (NCHS-1NC-1 )) 20,20, 

LCPERR * i 

WRITE (A. 9003) 

To continue 

IF (DD!2>> 30,30, 

READ <S. 1005) • YDD t 1 ) » 1 = 1 ,N X ) 

_ » R 1 TE Jl_6 • 2019) ( YP D< 1 1 *_I • N X ) _ 

~ 30 CONTINUE 

JEAD X ANO y COORDINATE LOCATIONS 
READ (5,1D05MXX(I).1*I,NX) 

READ (5 1 100SMTTI ,NY) 

WRITE <6,20051 (XXI I >. 1*1 , NX) 

WRITE « A. 2006) <YY< J) , J* 1 ,NY) 

READ THERMOCOUPLE CHANNEL NUMBER "AND 
CORRESPONDING X-Y COORDINATES 
READ (S, 1007) I IX, fr.NTCHl I X , I Y") , N* 1 .NTCLOC) 
WRITE (A, 2007) 

W R 1 T E ( 6 , 20 2 2) ( JEQ , J • J* l , N Y ) 

DO 90 1*1, NX 

WRITE ('6,2008) I , (NTCHl I , J) , - j*l ,NY) 

90 CONTJ NUE 

READ START STOP TIMES OR EVENT LEVELS ' 

IOPT * I O P Y 1 ♦ 1 
GO TO ( SO , 60) , I OPT 

.50 READ IS, 1003, END*A9n>TSTRT,TST0P. STATUE. NTME 
IFTnTSL.GT.O) WRITE I 6, 2021 ) 
write IA*2O03)TSTRT,TST0P,STATHE*NTME 
60 YO 70 ’ “ 

60 REAP 15, 1009, END-690)CESTRT,CEST0P,STATHE, NINE 
IF1NTSL.GT.05 WRItEIA, 2021 )' ~ 

write U,200«*)CESTRT»CEST0P,STATME,NTME 

70 IF 1 NTME«"l"E ,NLS) GO TO BO ~ 

! *_NLS_ _ 

WR I TE~1 A ,9002) I -- - - ~ 

LCDERR ■ 1 

bo continue " ' v. 

IF ( LCDERR) , , 73 D 

INITIALIZE AND DEFIJIE CONSTANTS 
C0N6 » HHO • CP /29. 

CON 12 * RHO • CP / 12* 

CON29 * BIGK • 29. 

- rp~AR — * u- - - — • 

NSCN * 0 

LTEMP * t 1 N Y*R ) / 10) • MX > S 
DO 90 1*1, NCH5 
DO 90 J*2,5 
90 TEMPI I , J) • n. 

IESTRf~S CESTRT ’ ' " 

I ESTOP * CESTOP 

npt ■ b 

LINES * 1000 
I YDD * NY ♦ 1 
IF (0012) > 110,110, 


3 



00424 

00425 
00430 

205* 
~ 206*" 
207* 

- — 

00 432 

208* 

100 

00434 

209* 

1 10 

00 4 35 

TlOi" 


00436 

211* 


00436 

212* 

C 

00434 

213* 

c 

00436 

214* 

C 

00441 

215* 

120 

00450 

216* 


0045! 

217* 

130 

00453 

218* 


00454 

214* 

140 

00457 

220* 


00460 

221 • 


00461 

222* 

150 

00462 

223* 


00463 

224* 

160 

00466 

225* 

170 

00467 

226* 


00472 

227* 

180 

00472 

228* 

C 

00475 

224* 

190 

005D0* 

i3o* 


00501 

231* 

200 

00503 

232* 


00504 

233* 


00505 

234* 


00507 

235* 


00510 

234* 


00511 

237* 


00514 

238* 

210 

00516 

239* 


005 1 7 

240* 

2 20 

00520 

241* 


00527 

242* 


00530 

243* 

230 

00531 

244* 

■ 

00533 

245* 


00536 

246* 

240 

00537 

247* 


O0S40 

248* 

250 

00543 

244* 

260 

00544 

250* 


00547 

25|* 

270 


252* 

253* 

~rsm — 

255* 

- 25 *. — ' 


00564 

257* 

300 

Oti^ra 

" 25¥*' - 


00570 

259* 


005T1 

' * 260*' ' 


00572 

261* 

310 





I ■ NT - 1 v 
DO JOO" j * ! 

1 F I Y Y I J ) . LE « YDP { t ! .AND,YY( J+l ) .ST* YDO( l ) 

CONTINUE ' ‘ 

_CO N T I N U F 

TIME (I ) * t I ME t 3 ) ; 

IFJNTSL-1 ) *150,150 

READ INPUT T APE _____ 

READ j LU I N , ERR» 1 3Q , END« 1 HQ 1 TIME! 1 1 » tTEMP 

CO TO 150 

WR ITE _I6*4O04) 

60 TO 120 ■ 

WRITE (6»4005)LUIN 
I .EOF ■ t 
60 TO HOP 

CONTINUE 

60 T O I 160*170) , I OP_T^ 

I F If I HIT 1 1-TSTRT l' 120,1 9 Q ,f 9 0 

IEVNT « TEMP ( NCC t 6 > ♦ 0*5 

IF I IEVNT- I I E5TRT + 2 )T” 1 80 , J 90*120 " 

IF ( lEVNT-t 1ESTRT-2) 1 120,190,190 
STARTING TIME FOUHO" 

DO 200 I«1 . NCHS 

Temp n;n » tempi i , 6 ) ■ "* 

TEMP U*2J “ TEMPI I, 6 I 
f IME I 2 7 * flME(l) ’ 

NSCN » j _ _ 

IF ( NTME • LE • 2 . OR • T I ME < t 1 ♦ LT . ST ATME ) GO TO 
NPT ■ NPT ♦ 1 

T I M I NPT 1 * T I N E (I J 1 “ 

DO 210 I»l 

TMP t T , NPTT; s^'T^riP I I , 6 > 

GO TO 330 

CONTINUE ' ’ - - * - - - - - 

READ lLUIN,ERR«23n.END-l*«0>TIMEl3) , (TEMPI 
GO TO" 290 ‘ * 

IPAR IPAR ♦ l 

WRlfE (6,40047 "" “ 

IF (IPAR-NPAR1 220,220,730 

CONTINUE 

GO TO (250,260) ,IOPT 

IF ( TIME 1 3 f-f STOP ) "2613,280,400 " 

IEVNT - TEMPINCC.6) ♦ .5 

if n ew-tteit op > 2 > t ztv v 270, 200 

IF ( TEVUT-I IEST0P-2) ) 200,400,100 

TIMET 21 « TIME ( 3 3 

NSCN ■ NSCN ♦ I 

IF <NSCMi27 300",, 300 ~ 

DO 290 I — I » NCHS 

TEMPT 1,3* * TEMPI ITS! — 

C0NTIN1IF 

DO 3 l 0 I * l , f! C H S ■ ■' 

TEMPI 1,4) * TEMPI I ,5* 

YEMpTT ,51 * TEMPI I ,6) 

TEMPI I ,1 ) * TEMPI!, I) ♦ TEMP11, 6) 



•s. 

00577 

00576 ' 

00577 

262* 

263* 

267* 



IF(NTME.LE.2.0R.T!ME( 
NPT * NPT + 1 
DO 320 I»1.NCHS 

0060 2 

265* 

320 

TMPC I , NPT 1 * TEMP(I,6 

00607 

266* 


TIM(NPT) * TIME12) 

00605 

267* 

330 

CONTINUE 

00606 

268* 


I F | NOTH 1 , ,365 

0061 1 

269* 


IF <LtNE5*LTEMP-55) 3 

00617 

270* 


WRITE (6.20001T1TLEI , 

00620 

271* 


LINES * 3 

00621 

272* 

370 

CONTINUE 

00622 

273* 


WRITE 1 6 • 20 1 3 1 T I ME I 2 1 

00625 

277* 


WR I TE ( 6 , 2023 ) (JEq.J. 

00637 

275* 


LINES * LINES ♦ LTEMP 

00635 

276* 


DO 360 1*1 .NX 

00470 

277* 


00 350 J*1 .Nr 

0047 3 

278* 


NN * NTCHl I , J) 

00677 

279* 

350 

WORD ( J J * TEMP ( NN f 6 ) 

00476 

280* 


WRITE (6,2010)1, (WORD 

00655 

28 1 • 

360 

CONTINUE 

00657 

282* 

365 

CONTINUE 

00657 

283* 

C 



GO TO 330 


.370, 


001207 
001222 
001225 
001235 
001237 
0012*12 
0012*12 
0012*1 *1 
001251 
00124** 
001247 
001247 
001303 
001323 
001335 


091335 

001335 


001334 

0013*12 

001343 

001343 


00457 
00657' 
00660^ 
0066 2 
00665 


28*1* 
285* 
28 6* 
287* 
288* 


(71 

I 

"4 


00666 

00672 

00477 

00700 

00700 

00700 


DO ERROR ANALYSIS IT »NT H£» POINTS HAVE B EEN SAVED 

IF(NTME»LE, 2«0R ,NPT»LT.NTHE) Go TO 220 

IE <LINE5+2*LTEMP-59> 370,370. 

LINES - 3 


001363 

001343 


287* 

270* 

271* 

272* 

273* 

277* 


WRITE (6.20001TITLE1 .TITLE2 
37 0 WRITE 1 6 » 20 1 5 ) NPT , T I N ( 1 1 iTir(NPT) 
LINES • LINES • LTEHP 3 

- CALL ERRFITlf 3951 


OUTPUT HE AT IN 6 R ATE - QQOT *Np PERCENT EllROR - ERR 


001363 

001363 

001700' 

001706 


001 R 1 Q 

00172*1 


00700 

295* . 

C 


~ 1 -' ' r ' • - * 

00701 

296* 



WRITE(6,2012) 

00703 

297* 



WRITE (6, 2023) '( jEq V J « J* 1 * Ny ) 

30712 

298* 



DO 380 1*1, NX 

00715 

299* 



WRITE (6,2010)1 , (QOOTI I »J> , j-t ,NY > 

00727 

300* 


380 

CONTINUE 

00726 

301* 



WRITE ("6.20 161 

00730 

302* 



WRITE(6,2023) ( JEQ , J , J* 1 » NY ) 

09737 

303* 



00 390 1*1. NX 

00772 

307* 



WRITE (6.20 10) I , (ERR< I , J > » J* 1 , NY ) 

00751 

305* 


390 

continue 

00753 

306* 


395 

CONTINUE 

007^7 

307* 



T!M( 1 > * T I M C NPT 1 

00755 

308* 



00 396 I ■ * 1 , NCHs 

057(50 

309* 


396 

tNPt T, 1 ) * THP 11 , NPT 1 ' 

00762 ■ 

310* 



NPT * 1 

00743 

ill* 



UO“TO "2"20 

.00763 

312* 

c 



00763 

5T3* 

c 


_ The — TThe Slice has been completed 

00763 

317* 

c 



00767 

315* 


700' 

continue 

■ 00765 

3 16* 



IF(NTME.LE.2.0R.NPT.LT.3) GO TO 77 

'00767 

317* 



IF (LINE5*2*LTEMP-59) 710,, 

■ 00772 

318* 



WRITE (4,2000)TITlE1 ,T1TLE2 

4 ... 

$■ t 



, 


- - 


001 * 13 * 

001772 

001772 

001772 


001772 
001775 
001752 
001505 
001505 
001527 
00 1527 
001531 
001552 
001552 
QO I 572 
001572 


OR Afl E~K D~OF FIX e' E N c 0 U N T E RETT " 


001572 
001606 
00 1 60$ 
ooj 6 1 0 
001612" 
001612 
□01612 
001612 
001617 
001617 
001631 
001637 


00776 

319* 


00777 

320* 


01001 

321* 


01001 

322* 

— c 

01001 

323* 

c 

01 001 

321* 

c 

01007 

325* 


01010 

326* 


01012 

327* 



01021 

01021 


01113 


01114 

01117 

OV01I20 


I 


00 01131 
01131 


01131 

01131 


01133 

01136 

01111 


“10 CONTINUE 

• Rj TE 16.201 5 ) NPT , 1 1 ITT iTt I liiW 
! F (NPT.LT, NTME) WRITE C & . 2 0 1 7 ) 


328* 

329* 


COMPUT E AND PRINT LEAST SQUARES HEATING RATE (QDOT) AND PERCENT 
ERROR (ERR) r “" “ ' - — ■ > — — 

CALL ERREIT(SIIO) 

«SR | TE ( 6 , 201 2 ) ~ 

• R 1 TE ( 6 , 2023 ) ( JEq . J , J* 1 .Ny) 


DO 120 I»l.NX J 

WRITE (6.2010) I » ( (jOOT ( I , J) , j» 1 , NY } 


01033 

01035 

330* 

331* 

120 

CONTINUE 

write (6.2016) 

01037 

01016 

332* 

333* 


WR I TE ( 6 . 2023 ) (JEQ.J. J*1*NY> 
DO H30 1*1, NX 

01051 

T31* 


WRITE (6.2010)1 , (ERR< 1 . J I . J* 1 , NY ) 

01060 

335* 

13 D 

CONTINUE ^ 

01060 

336* 

C 


01062 

337* 

110 

CONTINUE 

0106 3 

338* 


IE (NSCN-3) 680,, " 

01066 

' 339* 


NTSL - NT5L * 1 

01067 

315* 


LINES * 1000 

01070 

311* 


WRITE (6.2000)TITLE1 .TITLE2 

01071 

312* 


WRITE ( 6 • 20 1 8 j NSCN , T t ME M I . T 1 M Em 

01 101 

313* 


• RITE I 6, 2023) ( JEQ.J. J-l »NY) 


001 653 
001 6S 3 
001661 
001661 
001661 
001661 
001676 

ooTtoi 

001706 
001711 
001711 
00 1760 
001760 
□01765 
002006 


002006 

0 02026 

002026 

002026 

002026 

002031 


002015 

002017 


002061 

002071 


31S* 


316* 

317* 


DO 150 J*1 . MY 

NTCH( I , JT 


NN 

TEMPI NN , 1 1 


■ TEMP(Nn.I) / NSCN 


318* 150 WORD(J) * 'TEMP(NN.l) 

01122 319* WRITE (6.2010)1 , ( WORD ( J ) ,J»i ,NY> 


002122 

002122 

002122 

002123 


350* 

351* 


160 CONTINUE 


352* 

353* 


PROCESS LATERAL HEAT LOSS IN X-DJRECTION (EQUATION JT 


351* 

355* 

356* 


DO 170 1*1 .NX 
D0_170 J*UNY 
O X ( I . J ) « 0 « ~ 


01112 

357* 


QY ( 1 , J ) * 0. 

0TT13 

J51* 

170 

oboTf r» JT~ * o. 

01116 

359* 


IE (NX-2) 530.530, 

01151 

360* 


I DO * 1 

01152 

361* 


DO 510 .1*1, NY 

bus's 

362* 


1 GOOD » 0 

01156 

363* 


I F ( J.GT. IYDD) IDD > 

01160 

361* 


It * 0 

01161 

365* 

180 

11*11+1 

0IT62 

346* “ 


IE T I 1+2-NXl V*5 1*0 

01165 

367* 


NN I a NT CH ( 1 1 , J) 

01166 



ir (nni ) 180,1 so v~ 

.01171 

369* 


12-11 

ST 172“ 

3T3T* 

19F 

ir - * 12 + i 

01 173 

371* 


IE ( 1 2+ 1 -N X ) ,.510 

oira 

J72* 


NN2 * NTCH( 12, Jl~ 

01177 

373* 


IE ( NN2 ) 190,190, 

OI202 

371* 


13-12 

01203 

375* 

500 

13 * 13 ♦ l 


002127 
002131 
002155 
002155 
" 002154* 
002155 
002155 
002155 
□02155 
002155 
002156 
002161 
002170 
0021 75 
002175' 
002176 
002205 
002207 
00221 1 
_ 002215 
002222 
002221 
00i2f7~ 
002231 
' 002235 
002212 
002il1 
002,217 


D I 2 0 15_ 
01207 
012143 


376* 

377* 

378* 


1F_1 I3-NX1 , ,5 1 0_ 
N N 3 **■ NT CH{ l 3 , J ) 


00225 ! 


002258 


01323 


*•20* 


01328 
01327 
01332 
01335 
0 1 380 
1701381 

821* 
822* 
823* " 

828* 
825*" “ 
826* 

01381 

827* 

01381 

828* 

?01 382 

829* 

0X38 3 

830* 

rOi 388 

43 1* 

01388 

832* 


IG00D ■ IGOOD ♦ 1 
IF ( J2-ITDD J ,T80,S70 
IF (J3-1TDDS 575,575*580 
570 IF X j2-<TXDD*i ) )~59F»5907~ 

IF IJI-IYDD) 590,590, 

IDd "■ 2 ' ' 

575 QY1 I , J2) * C0N28 * DD(IDD) 


002527 

002532 

002536 


GO TO 600 
580 YYO * YDD(I) 

QYtI.J2) *= C 0 N 2 8 


___ w VY- JIH * < |TEMp(NN2,l » - 

TFIIPINMI ,rr> f ( YYTJ 2 T - YYUlFT - ( TEMp 1 NN 3 , 1 1 - ““ 

TEMPI NN2 , I 1 I / (YYIJ31 - YY(J2*)I 


OtJCD / (YY(J3)-vY(jl ) i • ( I TEMp I NN2 • 1 J ' » ' 
TEMP ( NN 1 , 1 ) ) / iyy ( J21-YYI J1 > ) - DP I 2 1 « (TEMP ( NN3 , l 1 - 


002588 
002550 
001558" 
002557 
002557. 
002557 
0026 12 
002618 
0026X5“ 
002615 


s 01213 

01218 
01218 

\ P L?!J_ 

379* 

380* 



I GOOD * IGOno ♦ 1 " — - — — — 

qx < 1 2 « J > « C0N28 • OD 1 IDD) / ( X X < 1 3 > - X X < I 1 > > * ( ( TEMP t NN2 • I > - 

UOctbt 

002268 

002267 



381 • 

382* 



• TEMP ( Nn 1 , 1 ) ) / ( XX ( j 2 1 - X X ( I 1 > ) - (TEMP(N N 3,1) - 

• TEMP(NN2j 1 )J/ ( X X_t 1 3 > - XXfl25»j 

002267 

flQ?267 



012)5 

01216 

383* 

388* 



QX t I 3 . J 1 * QX ( 12, J) 

!F( IGOOD.Eq* 1 > OXUliJ) * 0X112, J) 

002322 
□ 0237-3 



% 01220 
01223 

385* 

386* 



IF ( 1 3-N X ) ,510,510 - 

NN 1 * NM2 . ! 

002337 

002383 



v 01225_ 

38M 



NN2 ■ NN 3 . — — — — 

I 1 « 12 

002385 

002387 



01227 

390* 



12*13 - - 

GO TO 500 

002351 

002353 



% ul 2 Ju 

01232 
01238 
x 01283 

39! • 
392* 
;J*J« 
398* 


5 1 0 

CONTINUE 

WRITE 16*2009) 

002356 

002356 





RKITtl6| 2DZ3 ) l JEQ t J 9 J* 1 f N Y ) 
DO 520 j«l .NX 

002363 

002807 



01255 

396# 


520 

WRITE 16*2010) 1 , lox 1 I * J> * J*1 ,NY 1 ~ “ 

CONTINUE 

002809 

002827 



1 01*57 

01257 
01257 
x 01257 

01260 
01263 

397* 

398* 

399* 

800* 

901* 

902* 

C. 

53o 

CONTINUE "■■■— — 

002827 

002827 



C 


PR0CE55 LATERAL HEAT LOSS IN Y-D1RECT1 ON 

002827 

002827 





00 610 1*1 .MX 

* 002927 

002996 



\ u I 216 

, 01267 

903* 

909* 



IGOOD » 0 ' 

1 00 * 1 

002886 
□ □?U<17 



0,271 

01272 

01275 

905* 
_ 806* 
807* 

' : 


J 1 * 0 

ill * Jl ♦ 1 

002851 

002853 



908* 



Jr »Ji+2*Ni 1 910 *0 

NN 1 * NTCH ( I * J| 1 

002955 

0n29Al 



□ 1Z/6 

01301 

01302 

01303 
0L3Q6 

809* 

910* 

9 1 i * 
9) 2* 

~H j 3* ~ 



IF (NN1 ) 580,580* ; - — 

J2 * Jl 

002866 

002870 





550 

J2 ■ J2 * 1 _ 

IF (J2*1-NY) ,,610 * 

002873 

002875 



01307 

918* 



IF ( NN2 ) 550,550* 

002501 

002506 



01313 

81 5* 
816* 


_560 

J3 * JZ — 

J3 « J3 ♦ 1 

002510 

002513 



01317 
01320 " 

“ I ' • 

9 | 8* 
819* 

--- 


IF 1J3-NY) ,,610 ' — 

NN3 * NTCHI1.J3) 

if iNNa'rssa.sib". * — — — 

002515 

002520 


- 


5-10 


! ~ 


Oil^_ 

01314 

01345 

0134* 

01347 

01347 

01347 

01350 

31351 


01353 
01 354 


01357 

013*0 


013*1 

013*2 


013*3 

013*5 

on* r 

0137* 

orwoi 

01110 

01112 

01H2 


01112 

01112 

019)3 

01111 

61417 

01120 


01125 

0112* 

01127 

umr 

01132 

01133 

01134 

'oTinr 

01112 

umr 


11 3* 
'131* 
135* 


13** 

137* 

138* 

_439*^ 

110 * 

111* 

112 * 

113* 


590 


■ TEMP ( NN2, I 1 ) / ( DO t 2 1 • 1 YYD - YY(J2)) * OD ( 1 ) • ( Y Y < J3 > 

■ YYoTTT ‘ — ; 

GO TO 600 

yyd * yd d rr> * — 

0Y» I . J2) * C0N24 * OD 1 2 ) 


rrffpfNNr 


TH 

» YYOJJ _ (TEMP ( NN3 , 1 ) - 

*00 QY ( ( • J3 ) • (JY C I ,J2I ' ' 

1F1 IGOOO.EQ* 1 1 gY(t,Jl) « Q Y ( I . J2 ) 


/ CVYt J3>-YY( J1 ) ) • <DD( 1 >•( TEMP1NN2, 1 1 - 
f (DDT2 J *1yyO - Y Y ( J 1 » )~' ” 

TEMP(jgN2, 1 ) ) / 


DOTl J *TY Y ( J2 1 
( YY< J3)-YY( J2» > ) 


111* 

115* 


11 ** 

117* 


118* 

119* 

ISO* 

151* 

1ST* 

153* 

151* 

155* 


*10 


*20 


IF (J3-HY) 

NNI * NN2 
NN2 » NN3 
Jj ■ J2 
J2 - J3 
GO YO 5*0 
CONTINUE 
WRITE K.20I1I 
UWTf ErtT2023 r 
DO *20 1*1 .NX 

write mn 

CONTINUE 


1*10, *10 


< JEg » J • J-f.Hf) 

rrrjTij-T .nti 


630 CONTINUE 


15** 

157* 


~1SS* 

159* 

9*0* 

9*1* 


COMPUTE g (ARRAY ERR IS USED) IF W A NO Q Y NO T BOTH ZERO 

TTmeTTTT' * “INSCN-Ii / INSCN-2J 


~TKT~S 


zms- 

*50 I • 1 i 


~r rrrffEi 

NX 


TT 


9*3* 

~nv*~ 
965* 
9*** 
4*7* 
HI* 
4 (If* 


00 

“l ot) - 1 

DO *50 J • 1 * N Y 

nrrrjTsrvtYDoi ioo 

NN ■ NTCH ( I , J ) 
“ERR ( I , J) ■ a, 

IF (NN> 650, *50, 

Tinm 


TT 


470* 

471* 

“472* 


640 


.J) ■ UC * DO < IDO) 
TEMP ( NN » 3 I • 

roiitrr.jii *io,,aio 

*40,650*610 


T (TLHPIN N.bl * T eHHINN.4) - TERRI NN, 21 


01411 173* 

01997" 9749 — 

0115 1 475* 

’UT9475 4 mr 

014*3 477* 

01472 


*50 


ERR ( I , J ) ♦ QX ( j 
» gERTrHRTI .jl ♦ 


jl ♦ 9Y(I,J) 
TjUUTU.jT 


01474 
0T47* 
01505 
fffSIO 
.0 1^1 7 
0I52T 
01521 


~wi 

479* 


WO** 

481* 

4TY*~ 

483* 

~4lTlir 


6*0 


670 

TWO 


IF 

IF ( 9 Y ( ] iJH 
CONTINUE 
QD OT ( 1 tJ) « 

CONTINUE 

~write“ (AVTori) •— - 

WRITE!*, 2023) ( JEg , J • J* 1 »NY I 
ocr 6*0" in ,nx ~ — - 

WRITE ! 6 » 20 I 0 ) I , ( ERR < I * J > » J* 1 * NY ) 

CONTINUE 

WRITE (*.2012) 

*R I TE ( 6,2023) { jeq , j. JW1 . NY) 

00 670 1*1 , NX 

WW i TE *' (T720 1 OT I , ( gTToT f j , JT ji 1 , NTT* 
CONTINUE 


485* 

01574 Ha* c 

01524 487* C 

01521 ~9T«* ~C~* 

01525 989* 


TriTroFT- ; ,-pro — 

GO TO !50,*0),I0Pr 


00 2615 

002615 

002654 


002656 

002*57 

002*57“ 

002657 


002717 

002726 


002743 

002747 


□ 02751 
00 2753 

"002755 

□ 02757 
0027*1 

□ 027*2 
0027*7 
003017 
003Q 1 7~ 
003037 
00*037 
003037 
003037 
003037 
661037 
003072 
661672 
00307* 
00307* 
0031Q5 
OOlllS 


003111 
691 1 1 4 

003||4 

001124 

00312* 

603131 

003131 

661111 


COMPUTE and print meins of mean heating RATES* 
690 CONTINUE 


003171 

6612 1 0 

□03210 

6ol2l5 

003236 

00323* 

003256 

0032S? 

003260 

003T6TT 

0032*0 

661260 

003270 


oo 



003143 

—g M . 

_ _ 

□03150 




99o* 

99 r*~ 

992* 

993* 

_ 999* 
*195* ' 
_ 994* 
997* 
998* 
999* 
500* 
501* 
502* 
503* 
S09* 
SOS* 
5 04* 

*“ 507* ' 
SOB* 
509* 
SIS* 
Si !• 
S12* 
SI 3* 

S 1 9* 

5 15* 
S14* 

S 1 7* 
511* 

S 19* 
S20* 


write (6*200Q)TITlE 1 *TiTUE2 

WRITE ( 6 iTET2oTnT S L “ 

WRITEI4.2023) I JEq » J j J* 1 • N Y ) 

IF ( NTSL ) 730,730, ' 

DO 710 1*1 .NX 

D0“ 700 J* 1 . N Y ' ' ’ ~ 

700 Q DO THM ( I » J ) * qOOT WH ll.J l / NTSL _ 

WRITE (4.2010) I , (QDOTMM(T.J) ,J-UNY>~ 

710 CONTINUE 
720 CALL EXIT 

730 CALL~CHE*R “ ~ ' 

*» 

: FORMAT STATEMENTS 

1000 FORMAT ( 1 2A4 1 

1001 F0RMAT(SF10*3.213) 

1 002 FORMAT (7151 ~ 

1003 FORMA TI3F1 0.3, 15) 

tew format ( 2 f ro iorrr it. stk st - — 

1005 FORMAT (5F10.31 

1007 FORMAT (5(213,19)1 — 

1008 FORMAT ( l 5 X » A 5 ) 

2000 FORMAT .( 1H1 / 14X , 1 2A47I4X . 12A4) 

2001 FORMAT! I H O , 1 S X . * MASS DENSITY OF TEST MATERIAL (RHO) 

' , 5 X . 1 SP E CTFI C HEAT OF TES T MATER fAL ( CAT 


01409 

522* 

01409 


01609 

529* 

01409 

52S* 

01409 

526* 

0T4TJ9 

527* 

0140S 

528* 

~Ol405 

5J9* 

01405 

S30* 


0TSI5S — 
01 604 
ortti4" 
01404 


-5TT* 

532* 

533* 

539* 


1 4x » * thermal conductivity (bicki 

. SX . • TH I CFNESS of TEST MATERIAL (WIT 
14X . ’thickness of test material iddzi 

' • SX. 'NUMBER OF X - C OoSTTTN ATES ( HTT 

1 4X i ' NUMBER OF Y - CO OR 0 INA TES (NY) 

14X. 'NUMBER OF CHANNELS ON TAPE (NCHS) 
.Sx.'commAno CHANNEL NUMBER ( NCt ) 

1 4X . » T I ME INTERVAL OR EVENT CODE ( I0PT1 ) 

» 5X » • ACCEPT ABLE NO • OF PARITIES HlpAf ) 

1 4X . ' T I ME HISTORY TAB OPTION (NOTH) 


FORMAT! 


14X. 'NUMBER OF CHANNELS ON TAPE 

75x . 'command channel number 


(NPAR) 

(NOTH) 


FORMAT ( |H0, 15X. 'START t!*E •'.F10.3,' 


STOP TIME ■' 


.'LEAS T SQU ARES COMPUTATIONS ( N THE ) *»,Ho) 
20 0 *) TORN ATTTWDT ISXT'To H^N6~ETFNT "STtOTT LEvELT 

. . ' *» . F I 2 . 9 / 

• T4x.' r LEA5T SQUARES START TIME ■> T «FTO. 3,' 

.'LEAST SQUARES COMPUTATIONS (NtME) *',Il0) 


>• ,F12.9 
■' . F 1 2. 9 
»• .F12.9 
rr,F 12.9 
■' »F 12.9 
>'.112 
»« .112 

•'.112 

■t;tt2 

■'.112 
■• . 112 
■',112 

1 . F 1 0 • 3 / 

ITTfS For ' 


01412 
- 014 !X 
i 2 01419 

~rorsT5- 

1(01414 
01617 
s 0 1 4 2 0 
7 01621 
601421 


592* 
59T* 

599* 

* 595* 
594* 


200 4 FORMAT (1H0 .1SX , • Y-COORD I N ATES . . . • ( /2 I X . 1 OF 1 0 , 3 )) 

■ ZDOT FORMAT ( I HOI TS T, ' I NPU7 TTPE CH)TWNEt ‘ NOmFTE R s'“.“ » T T )~ 
2008 FORMAT (IH , 1 5 X »• I • 12.101 10. (/22X, 101 10) ) ‘ 

2009“ FORHTTTTH0TIOX EKN t ATERAL HERTHLOSS TN “X 1? OTRFrCr; 
2010 FORMAT (1H . 5 X . ♦ I * « . 1 3 . I OF 1 0 » 3/ ( 1 2X . 1 OF 1 0 . 3/ ) ) 


2012 F ORMAT ( 1 HO. I OX . 'MEAN HEATING RATE 

2013 FORMAT TTHCT. TUX , 'TIME- t 7fTo* 3/7I I XT' £0 R R EttT“TEHT> E« A f UR E . ..1 
20)9 FORMATdHO.lOXr'MEA'N HEAT FLUX WITHOUT CONDUCTION CORRECTION .. 
2015 FORMAT TTHOTI OX, 'LINEAR CURVE" FTT USING - '» fTf POI’NTS.' - TTME - " i 

.FlO. 3,* TO*, F10. 3,* SEC') 


4 

3 


01622 

ou«~ 

01 629 

oiAyr - 

01625 
01625' 

01626 
01626 
01627 

01630 
0163 0 

01630 

01631 
01631* 

01631 

01632 

01633 
oiiW 

01635 

01636 * 
01636 

' 01636 
01636 
01636 

01636 
0163T 

01637 


rvj S77* 

Q I 61.2 571* 


01615 SCI* 


2016 format ( 1H0» 10X, ‘PERCENT standard ERROR ..,*) 

701 7"FBfiW*T — TTHr iTBTyrs** FEWER~rOrNT5~Tff6N ~9E Q u E5TETB '** • — 

2018 H °» *o*» * I E " P 5”* T ^ RE:S averaged oveR* t i 5,* points* time. in» 

2 0 19 FORMAT ! 1 HO , 1 S X . * Y-COORp INATE5 OF THICKNESS DIFFERENCE JUNCTION...* 

. »T7TTx 1 1 of lo, JH ““ " 

2020 FORM AT! 1 HO i 1 OX . « MEANS OF MEAN HEATING RATES F0R»tl3»* TIME PERIOD* 

«i*S...*> * ’ ~ ' 

2021 FORmATC l M 1 1 

2022- FORMAT : 

.( 28X.A3,I2,5X.A3,I2.5X,A3,I2.5X.A3,12,5X.A3,I2,SX,A3,I2.5X,A3.I2, 
.5X.A3.I2.5X.A3.I2.5X.A3, 12) r ~ 

2023 FORMAT 

.( 17X ,A3iI2,5X , A3, 12.5x763*. 1 2 . 5 X , A 3 , 1 2 , 5 X , A J . I 2 , 5 X , a 3 , J 2 . 5 X , A 3 . I 2^“ 
.5X , A3 * 12 .5X i A3 , IZ.5X • A3 »12> 

1001 FORMAT ( | X . * • •• ERROR - NX OR NY IS BAD*) 

1002 FORMAT 11 X . *••• ERROR - NTME IS LARGER THAN *tl3) 

1003 FORMAT! 1 X,'**» ERROR TCTQ M aNy~CHANNELS *1 

1001 FORMAT ( IHO , 1 OX i •• *' • PARITY ERROR RECORD DELETED • • •») 
"fTTOrTeRMATTYTOi lor.TTND Of TTLe '■.EWgT»fffTEffE6~gy"D TfrT~r"i"l 2 » ♦ PROCESS INS* 
..* OF ACCUMULATED DATA RILL BE ATTEMPTED.*) 






■•••••••••••••••••••••••••a 



01653 
in 653 

01655 
0T656T 

01656 

01657 
01662 
01X65 
01666 
0 1 670 * 


587* 
~5ff*v " 

S87« 

Wr’ 

591 • 

59 J* 

593* 

599*~ ' 

595® 

5966— 


INTERVAL OF DATA TO A FIRST DEGREE POLYNOMIAL AND TO COMPUTE 


1 SUM 1*0 * 0 


£ 3 n — i t in r | 

SUM X,( I I 

2 m sum 2 «r“TTMTFn 

SUM (XII) ••21 


125 CONTINUE 

fslJM - r(T) Tm42 - 
SUM6 « SUM29 • 2 

‘ 5UHS~*l!PT“i SDMI SUM 6 

50“" 2 oo“ r* ttnx — — 

DO 200 J = 1 .NY 
i Ob' *“1 ~ 

|1F( J.GT. I Y DO) IDD ■ 2 


01671 577* 


IF (NN.LE.O) 

GO T 

■ 01 673 57§* 


5UM 1 *“01X3 


1 01*79 577* 


SUM3 ■ 0.0 


01*75' “ 5rrTT» 


5055'*’ 0.0 


01676 601* 


DO 150 m * 1 

,NPT 

0TS7S 607* ' 

C 

SOM Yt I ) • Y ( I) 


SUM 1 + TIM(M) • TMPtNN.MS 


003907 

003107 

003107 

' 00 390? 
003907 

'oW 

003907 

QQ39Cf7“ 

003907 

00$907 

003907 

003907 

003907 

003907 

003907 

003907 

B03907 

003907 ' 

00 3907 

00 39! 7“ 

003907 

003907 

0 03907 

003907 

003907 

003907 

0039Q7 


003907 


907 
003907 


00390 

003907 


00391 

003916 


00399 

003996 

00J69S 

003996 


003956 
*003956 
003956 
003960 
003960 
003969 
003513 
003513' 
003515 
00 3525* 
003526 
003532 
003539 
003535" 
003595 
063595 
003595 




01701 

605* 

C _ 

SUM 

01702 

60S*" “ 


SUM 3 » 

01702 

606* 

C^ 

SUfl 

01703 

607* 


SUMS «' 


01.705 
0170** 
o 17 aj, 

01706 

0 1707 
01707 
0 1 7 l 0 _ 

01710 
01710. 
0171 1 

0 1711 

01712 
0 ,7 L* 
01713' 
0I7J5_ 
01715 
0171 6 
01717 
01722 

01722 

01723 
01725" 


__ 150 CONTINUF 

C ~ SUM xTf ) • SUM Y'm “ ' ' - - - — - - 

SUM 7 «» SUM 2 • SUM3 

C ' SLOPS ' : ‘ " ~ •*"_ 

SLOPE ■ <NPt * SUM1 - SUM7> / SUMS 

c sigma " ' ; 

SIGMA ■ 1 .0 / NP7 * SQRTI^NPT • SUMS - JUm3»*2 

1 - INPf • SUM 1 - SUM7 ) **2" 7 SUHliT " — 

C_ _ „ S.TAN0AP0 j:RR0R_ 

SERR * NPT * SIGMA / SQRTI <NPf-2> • SUM#) 

C PERCENT ERROR 

ERR ( I » J ) * SERR / SLOPE * 100*0 
_C HEATING RATE 

IJDOT • I » J 1 » CON 12 * ODflOD) • SLOPE 

GO TO 200 * 

175 ERRTirJl ■ 0.0 ‘ ‘ “ ” — — 

QDOTflsJI * 0.0 

200 CONTINUE . ~ : 

RJETURM. 

c — 

759 . BfRir - -i&i?999) 

9999 rORMAf (THO , 10x , •TOO EE* POINTS AVAILABLE FOR CURVE PIT*') 


01726 

630 * 

RETURN 1 



01726 

631 * 

C 



01727 

* 32 * 

END 

. . 



_ 003 S 5 S 
‘ 003552 
00 3 S 52 
" 00355 S 
00354 ^ 
003562 
003562 
003542 
003545 
003545 
003573 
003573 
003573 
053412 
003412 
0034 1 i 
00341 * 
003421 
003524 
003430 
003530 
003434 
003636 
003434 
003652 
003456 
003656 
003456 
003714 


END OF compilation;" 


i'agnoSticS. 


SCOPOUT" TPPS, ,tmns<n.“ ' 
PURPUR 0026 - 12 / 20 - 03:02 

‘TTPs r itc 

bprep 




V, 

7 


ATTACHMENT 6 


SAMPLE INPUT 






(TO BE SUPPLIED LATER) 
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ATTACHMENT 7 


SAMPLE OUTPUT 


(TO BE SUPPLIED LATER) 
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ATTACHMENT 8 


CORRESPONDENCE 



THIS FORM MUST BE COMPLETED BY TYPEWRITER 


01 4 

JSC 

ot 7 PROGRAM NO 

Q614 

JSC COMPUTER PROGRAM ABSTRACT 

01 14 DATE (MMDDYY) 

1-21-75 

01 20 TITLE OF PROGRAM (62 CHARACTERS MAXIMUM) 

SURFACE HEATING RATE OF THIN 
SKIN MODELS 

0 172 SYMBOLIC NAME 
(9 CHARACTERS 
MAXIMUM) 

THNSKN 

PARENT PROGRAM j 

02 14 

CATEGORY 

02 16 SITE 

02 19 

program no 


02 26 
CAT- 
EGORY 


02 27 
LANGUAGE 1 
NO* 1 

FOR 

V 


02 32 
LANGUAGE 
NO. 2 


02 37 KEY WORDS (8 MAXIMUM SEPARATED BY COMMAS) 

Arc Jet, Heat Rate, Thin Skin, Thermocouples 


WHOM TO C 

ONTACT ABOUT THE PROGRAM 

05 48 STATUS 

05 49 

A. THIS PROGRAM 
PI IS NOT FOR 

SHARING 
__ B. limited 
LJ SHARING ISEE 

ABSTRACT) 

OS 14 CONTACT (LAST 
N AME) 

Stephenson 

05 28SITE 

JSC 

05 310RGN 
CODE 

FD57 

05 39 PROJECT NO 

4300 

05 45 NASA 
CENTER 

0A. UNDER DEVELOPMENT 
Q B. OPERATIONAL 
0 C. COMPLETED 


DATES 


05 58 REVISION CODE 


TIME AND COST FOR DEVELPOMENT 


0 5 50 

INI Tl ATED 
MMY Y 

11-74 


05 54 

COMPLETED 
MMY Y 

12-74 


05 59 MAN-MONTHS 


Q A. REVISION 
0 B. CANCELLATION 


J L. 


.8 


59 60 61 62 63 


05 6 4 MACHINE HOURS 

I 1 2 i * I- 


64 65 66 67 68 


05 69 CO^UTER 
TYPE 

1108 


05 74TOT AL COST 
(DOLLARS} 


-L 


1 J 1. 


74 75 76 77 78 79 BO 


ELITE MARGIN 


CARD 

NUMBER 


3 ABSTRACT 

j 

THNSKN computes the mean' heat ing rate at up to 1 0 d 1 


locations on the surface of thin skin transient 


PICA MARGIN 


heating rate models. Output is printed in tabular 


form and consists of time history tabulation of 


temperatures, average temperatures, heat loss 

without conduction correction, mean heating rate, 
least squares heating rate, and the percent standard 
error of the least squares heating rates. The input 

. 1 • 1 1 1 . 1 Ti r t m n a »» 


error of the least squares heating rates. 

tape used is produced by the program EHTSQ3 
( SCR a3-15-10/.McEryd£.) 


RELATED DOCUMENTATION (66 CHARACTERS MAXIMUM SEPARATE EACH REF BY COMMAS} 


42 


ICP Fnrm Id? 1 Luo 71V 
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NASA- KC 



OPTIONAL. FOAM NO. to 
MAT ISC IsmON 

mmA ppm* (4t cm) toi-iM 


UNITED STATES GOVERNMENT 

Memorandum 


Cr— TnSsSSoa nnd Ansslyta D:vie:on 
fcil£A Marked Spzc^iJt Carter 
Kciistaa,. iszaa 7703a 


TO 


date: October 27, 1969 


: LEC/Manager, Test Data Reduction 
Department, 672-10 

In reply refer to: 

FROM . : ED5/Acting Chief, Data Processing Branch 11-69-622 


SUBJECT: Requirements for a Thin Skin Transient Heating Rate Model 

Program 


Programing support is required to develop a program to compute 
heating rates of thin skin transient heating rate models. 

This program will be developed for the Experimental Heat Trans 
fer Section of the Structures and Mechanics Division. The 
program will be used to support definition of the aerothermo- 
dynamic design of the MSC space shuttle vehicle. 


It is necessary that the program, specified by 
requirements document, be operation by October 
liminary requirements have been reviewed by Dr, 
project engineer, and forwarded to Mrs. Yvonne 
Work on this task should be charged -to project 
For further information, 
extension 5533. 


the attached 
31, 1969. Pre- 
Carl D. Scott, 
Chempinski . 
number 4556. 



Enclosure 
cc : 

ES56/Dr. C. D. Scott 


ED5 7 :E0Grice(LEC : JMRandal) :ddc 


10-27-69 


DATA PROCESSING REQUIREMENTS 
for a 

THIN SKIN TRANSIENT HEATING RATE MODEL PROGRAM 



Prepared by: 


' iJ/ /fatty, t 

.J.^.-ftadal TfiS&Th't 

(EEC Requirements & Analysis 



Approved by: 





K. J^Xaj^enJiferg J/O Z4 
LEG Requirements & Analysis 


Concurrence by: 


Concurrence by: 


e, o. a-w. i 


E. 0. Grice /c /:? -5/f/ 

CAD Requirements & Analysis 


A t 


n 




.r' < 






Carl D. Scott / .;/ pu - ■ 

Experimental Heat Transfer Section 
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PROJECT DISCUSSION 


The Experimental Heat Transfer Section of the Structures and 
' Mechanics Division has requirements for a digital computer 
program to compute heating rates of thin skin transient heat- 
ing rate models. These models will be tested in either the 
1.5 Megawatt or the 10 Megawatt Arc-Heater facility. The tests 
in which such models will be used will be directly related to 
current MSC space shuttle vehicle activities. The program out- 
puts will be used to support definition of the abrothermodynamic 
design of the MSC space shuttle vehicle configuration. 

II TEST CONFIGURATION 

The thin skin transient heating rate models will be constructed 
with a thin metal skin. Thermocouple instrumentation will be 
attached to the back side of the skin in rectangular arrays and 
used to measure temperature. When these models are subjected 
to a non-uniform constant heat flux by the arc-heater equipment, 
the surface temperature of the skin will increase according to 
the laws of heat conduction and capacitance. The skin will be 

t 

sufficiently thin so that the temperature gradient normal to 
the skin will be negligible. Therefore, it will be assumed 
that the thermocouples attached to the back surface actually 
measure the temperature of the skin at each location. 

The analog signals output by the thermocouple array will be 
digitized by use of either the 160 channel data acquisition 
system located in Building 222 or the 50 channel data acquisition 
located in Building 262. The digital data tape output by either 
, of these systems will be compatible with the CAD UNIVAC 1106 digital 
computer system. The existing EHTS03 (Q566) program will be used 
to calibrate the experimental data in engineering units and out- 
put a data tape for input to the required thin skin transient 
heating rate model program. See Figure 1 for a diagram of the 
data flow and test configuration. " 
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1.5 MEGAWATT 
ABC HEATER 


50 CHANNEL DATA 
ACQUISITION 
SYSTEM 


160 CHANNEL DATA 
ACQUISITION 
SYSTEM 


^mgitalN 


TAPE 


DIGITAL 


TAPE 


EHTS03 PROGRAM 




UNIVAC 1108 
DIGITAL COMPUTER 


'DIGITAL 
i TAPE 


J DATA | 

! TABULATION 1 


THIN SKIN TRANSIT 
HEATING RATE - 
MODEL PROGRAM 


UNIVAC 1108 
DIGITAL COMPUTE?. 


DATA 

TABULA! 


Figure 1 

TEST CONFIGURATION 









Ill 


SATA PROCESSING REQUIREMENTS 


' The hardware required to accomplish the data processing task 
Will be the UNIVAC 1108 digital computer. The software re- 
quired will be the existing EHTS03 (Q566) program and the 
Thin Skin Transient- Heating Rate Model Program specified in 
Appendix A. 

The' data processing procedure will consist of the steps listed 
below: 

1 • The project engineer will have the data tape and the 
lead card input for the EHT 0 O 3 program delivered to 
the test data coordinator in Room 228 of Building 12. 

2. The data coordinator will issue instructions for 
computer processing and notify the project engineer 
when the processing is complete. 

3. Ihe project engineer will review the EHTS03 data and 
provide lead oard setups for the thin skin transient 
heating rate model program. 

• * * \ 

4. The data coordinator will issue instructions for 
computer processing and notify the project engineer 
when the processing is complete. 

IV . TECHNICAL REFERENCES 

1. Memo: ES 5/9-24(9) /249M, Memorandum to: ED/Chsef, Computation 

and Analysis Division; from: ES/Chief, Structures and Mechanics 
Division; Subject: Requirements for a computer program to 

calculate heating rates using thin skin models. 

2, Program Documentation, Experimental Heat Transfer Data Reduction 
Package #3 (EHTS03) , Program - Q566, Project 4556, August 1969. 
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APPENDIX A 


PROGRAMING SPECIFICATION 
for 

THIN SKIN TRANSIENT HEATING RATE MODEL PROGRAM 

I PURPOSE 

The purpose of this program will be to compute the heating rate 
at .each of a number of locations on the surface of thin skin 
transient heating rate models. Experimental data, calibrated in 
engineering units, will be obtained from the ABEL tape output by 
the EHTS03 computer program. The program output will consist of 
time history tabulations of computed data with options for input 
and intermediate computation tabulations . 

II INPUT 


The program inputs will be the ABEL time history data tape output 
by the EHTS03 program and the lead card inputs necessary to specify 
computational constants, input data channels, command channels, op- 
tions, and command event levels. The magnetic input data tape will 
consist of FORTRAN written data records in a single data file. The 
first data word in each record will be the time word.' The remaining 
words will be the data channel words ordered by channel number. The 
tape will be written with 800 bit per inch density and odd parity using 
the UNIVAC 1108 computer. 

The lead card inputs will include, but not be limited to, the 
following items. The exact card formats will be determined by the 
programer. 

1. Two cards, each containing at least sixty column alphanumeric 
entries will be used to identify the data contained in the 
output tabulation.. This input will appear in the title 
of each page of the tabulation. 
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2. Coordinate positions for each thermocouple with a maximum 
array of 50x50 positions. The upper left hand position 
will be designated T = T 1 1 as shown by Table 1 in 

, AJ I I 

Section III. 

3. Input tape data channel number for each thermocouple in 
T xy array * 

4. Physical locations of each row and column. All 
locations will be determined from the rc-w and column 
containing T^ = T^ in inches. 

5. Input tape data channel number for the command channel. 

6. Command channel signal level used to start computations. 
Millivolt level. 

7. Command channel signal level used to terminate computations 
Millivolt level, 

8. Input tape data channel number for the time channel. 
Normally channel 1 . 

✓ 

9. Start computation time. 

10. Stop computation time. 

11. Coordinate positions of known erroneous thermocouple data. 

12. Option to tabulate tape input temperatures. 

13. Option to tabulate intermediate computation results. 
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u. 


Computational constants : 

a. Hass density ( Qf .) of the test material in pounds 
mass per cubic foot. (ibm/ft^) 

b. Specific heat (C^) of the test material in Btu/d.bm °F 

c. Material thickness (d) in inches 

d. Thermal conductivity (K) in Btu/ft sec°F 


TABLES 


The thermocouple instrumentation used to measure the temperatures 
on the thin skin transient heating model will be arrayed in the 
rectangular coordinate system shown by Table 1 . The maximum size of 
the array- will be 50x50. The minimum array will be 3x1 or 1x3. 

Only line, rectangular or square arrays will be defined. Other array 
configurations may be computed if dummy data is input to fill out the 
array. Even though the thermocouple array will generally be rectangular, 
the spacing between rows and columns will not normally be the same. 
Therefore, the distance of each row and column from T 11 will be a 
separate card input. The array will be computed as for a planar 
surface regardless of the actual configuration of the model. That 
is, the last row will not be adjacent to the first row as in the 
case of a olyinderical surface, etc. Spherical and conical surfaces 
will not be defined other than as rectangular arrays. , - 



TABLE 1 

THERMOCOUPLE ARRAY 


8-9 


IV 


COMPUTATIONS 


Equation 1 will be the energy balance equation used to describe 
the heating applied to the model surface at each coordinate 
( X i , ^ t j) ^ or eac k < ^ a ^ a ^ me sample 

Equation 1 : 



where: 



q^..: Heating rate to the model surface 

P : Mass density of surface material 
d: Material thickness of surface material 
0^: Specific heat of surface material 



Derivative computed by the use of equation 2 

Derivative computed by the use of equation, 4 
Derivative computed by the use of equation 5 
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Equation 2 will be used to compute the partial derivative of 
temperature with respect to time for the thermal mass around 
the thermocouple test point for each data time sample (t^.) . 

Equation 2: 


a* 


U 


a t 


= \ + zo k \ 


•Where: 

b^: Determined by solution of equations 3a, 3b, and 3c 

C^: Determined by solution of equations 3a, 3b and 3c 

Equations 3a, 3b and 3c will be solved for b^ and by the use 
of the Gauss-Jordan or other appropriate method for the solution 
of simultaneous equations. 


T u <W = a k + Vk-i +c k ( Vi )2 

i 13 <v> = -k + Vk + c k ( v 2 


: T iJ^ t k+1 


> =a k + Vk+i + c i 


At-1^ 


Where: 

W°k : 
Vi : 

A 

Vi 1 


Coefficients of a second degree polynomial 

Time of preceeding data sample from ABEL input tape 

Time of current data sample 

Time of succeeding data sample 
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Equation 4 will be used to compute the lateral heat loss in the 
X direction for input to equation 1 for each time sample ( t^) . 



k : Thermal conductivity (card input) 

d : Material thickness (card input) 

X^ + .| - X^ : Distance between thermocouples located on either 

side of X^ thermocouple, etc. 



Input data temperature difference between thermocouples. 


Equation 5 will be used in the same manner as equation 4 to determine 
the lateral heat loss in the Y direction for each data time sample. 


Equation 5: 


•Y 

q£ .j = 24 kd 


3+1 


“ y i -1 


FT _ T 

1.1 i.J-1 

y j ' y J -1 


+ 


T - T 

-ij+1 - i ll 


3+1 - 


- 


1 

Equation 1 will be set to zero for the two end points in a test array 
consisting of a single line of thermocouples. Further, equation 1 
will be set to zero for the four corner transducers in a square or 
rectangular array. That is: 

2< i< 1-1 , 

2 < j < J-1 

where: I and J are defined as the greatest i and j limits, respectively, 

of the transducer array. 

Equation 1 will be computed for the test points in the transducer array 
that lie on the four sides of 3x2, 2x3 or larger arrays. Either 
equation 4 or equation 5 will be modified as required to prevent 
erroneous results due to lack of valid inputs for the transducers 
located on the sides of the test array. For test points located on 
the upper side of the test point array, equation 6 will be used 
instead of equation 4. 

Equation 6: 

For test points located in the lower 3ide of the array, equation 7 
will be used in lieu of equation 4* 
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Equation 7s 


i 

I' 



Equation 5 will be valid regardless cf whether equations 6 or 7 
are used for upper and lower limit test points. 

For test points located on the left side of the array shown in 
Table 1, equation 8 will be used in lieu of equation 5. 

Equation 8: 



For test points located on the right side of the array, equation 9 
will be used in lieu of equation 5. 

Equation 9: 



Equation 4 will be valid regardless of whether equations 8 or 9 
are used for left or right test points. Note that the end points 
of a linear array and corner points of a rectangular array will not 
be computed, as previously stated. 
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Card input will be used to specify which if any, of the input test 
points have invalid tape data input due to hardware malfunction in 
the test hardware or the data acquisition system. Computations for 
the test points on the four sides of the invalid points will be made 
by using the input temperature from and the distance to the thermo- 
couple on the opposite side of the faulty point. For example, if 
T 2 p were faulty (see Table 1), q^ ^ would be computed using Tp ^ 
for the temperature input and for distance input to 

equation 5. Note that the program will be designed to salvage data 
only for points within the array. A bad point at any edge or corner 
of the array will be eliminated by deleting that row or column from 
the program input. Substitution of data from a similar point may 
also be done by card input. No provision will be made for computing 
data around specified bad input points that are separated by at least 
one row or column of good input data. 

Computations will be controlled by card input star - ..-stop times or 
by card input command channel levels. Either one or the other may 
be used, but not both. Star top times will be input inseconds and 
milliseconds. Command Channel start and stop levels will be input at 
the same levels used in the ETHS03 program. These levels will be 
30,35,40,45,50,55,60,65,70,75,80,85,90 and 95 millivolts plus or minus 
2.5 millivolts. 
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PROGRAM OUTPUT 

* 

The program output will consist of tabulations of the input lead 
cards, thermocouple array format, computed and input data. No 
data tape will be output. Each page of the tabulated output will 
be identified by the card input identification title ii* addition 
. to any other identification, column headings, etc.,' The tabulations 
will include: 

1 . Tabulation of all input card data. 

2. Row by row tabulation of input tape data channels: This 

tabulation will start with the channel number for T = 
and list the input data channels in 10 columns across the 
page. See Table 1 for a typical format. Channel numbers 
assigned will be set to zero for those channels specified 
as erroneous by card input. 

3. Row by row tabulation of computed heating rate (q. .) for 
each time history slice. The time of each sample will 
precede each beginning row of computed data. Time will 
be output in seconds and milliseconds. 

4. Optional row by row tabulation of tape input' data tempera- 
tures for each time history slice as tabulated in 3, above. 

5« Optional row by row tabulation of computed X and Y heat 
losses (q^ and q^) as in 3, above. 

Note that tabulated data outputs for non-computed end points, corner 
points and bad data points will set to aero. 
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aeusTiN uitiruE ititeu 
16811 EL CAMINO REAL 
HOUSTON, TEXAS 77058 
TflEPHONt, (AREA CODE 713) 4U-00«0 

ELECTRONICS COMPANY 

In reply refer to: 
TDR#2247 


October 31, 1969 


TO: ED5/Chief, Data Processing Branch 

PROM: LEC/Manager, Test Data Reduction Department 

SUBJECT: Technical Memorandum 11-69-622 s October 27, 1969, 

Requirements for a Thin Skin Transient Heating 
Rate Model Program 



$ ^LOCKHEED 

v - . a' 


Subject Technical Memorandum was received by the Test Data Reduction 
Department on October 29, 1969 . 

Responsibility for the requested programming support has been assigned 
to Mrs. Yvonne Chempinski. Mrs. Chempinski has been working from pre- 
liminary requirements documents. The program has been completed on 
schedule with test output delivered to Doctor Scott. Two minor addi- 
tions were requested by Doctor Scott and will be completed shortly. 

Please direct any additional communications regarding this project to 
Mrs. Chempinski at HU 8 -OO 8 O, extension 345 or Mr. Len Martin, ext. 336. 

Unless otherwise directed, all effort expended on this task will 
continue to be charged to Project Number 4556. 



DLJ : JDG : pm 
Distribution: 


J. Fisher/ED57 
E. Grice/ED57 
J. Randall/Cl8 



L. Martin/Bll 


Pi Garrett/Bll 
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to : 

FROM : 

subject: 



lots- to. 


Computation ind A*aly§!t DitWo* 
f| ASA Manned Spacecraft Cantor 
Hraatno, Tuu 77058 

LEC/Manager, Test Data Reduction DATE: November 18* 1969 

Department, 672-10 

In reply refer to: 

ED5/Chief, Data Processing Branch 11-69-634 


Requirements for modification of the THNSKN (Q614) program 


OPTIONAL FORM NO. It 

MAY \m EDITION 

•SA r* MR (41 cm) |tt*1t.t 


UNITED STATES GOVERNMENT 

Memorandum 


Programing support is required to modify the Thin Skin 
Transient Heating Rate Model Program. The modifications 
requested by Dr. C. D. Scott of the Experimental Heat Trans- 
fer Section of the Structures and Mechanics Division are 
listed belcw: 

I 

a. Change the card input format for X and Y coordi- 
nates from F 5.1 to F 10.3.. . 

b. Change the data. output print format to F 10.3. 

c. Make the tabulation of tape input data temperatures, 
X and Y lateral heat losses as standard output rather than 
optional . 

d. Compute end point and corner point data for the 
thermocouple array rather than setting this data to zero. 

It is understood that Mrs. Yvonne Chempinski has started 
work on these modifications in response to a verbal request 
from Dr. Scott. Work on this task should be a continuing 
charge to project number 4556. For further information, 
please contact E. 0. 


cc : 

ES56/Dr. C. D. Scott 
ED57/J. L. Fisher 

ED57 rEOGrice (LEC : JMRandal) :ddc 11-18-69 
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LOCKHEED 


ELECTRON! C S 


MUST 

16 011 EL CAMINO REAL 
HOUSTON, TEXAS 77060 
TUCPHONCt (ARIA CODE 711) 418-0019 

COMPANY 

In reply refer to: 
TDR02267 


November 25, 1969 


TO: ED5/Chief, Data Processing Branch 

FROM: LEC/Manager, Test Data Reduction Department 

SUBJECT: Technical Memorandum 11-69-634, November 18, 1 969 

Requirements for Modification of the THNSKN (Q6l4) 
Program 


Subject Technical Memorandum was received by the Test Data Reduction 
Department on November 21, 1969 . 

Responsibility for maintenance of the THNSKN program has been 
assigned to Mrs. Yvonne Chemplnski. The program modifications 
specified in the subject Technical Memorandum have been completed, 
and the modified program has been validated. 

Please direct any further communications regarding this task to 
Mrs. Chemplnski, HU8-0080, extension 345. - " ; 

As directed, all effort expended on this task will continue to be 
charged to Project Number 4556. 


DLJ : JDG:pm 

Distribution: J. L. Fisher/ED57 

E. 0. Grice/ED57 
J. M. Randal/Cl8 
Y. Chempjnski/Bll 
J. D. Garrett/Bll 



s 

ION OF LOCKHEED AIRCRAFT CORPORA T ! C ,'i 


A 0 j V I S 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
MANNED SPACECRAFT CENTER 
Houston, Texas 77058 


IN REPLY REFER TO: 


11 - 69-636 


2 4 NOV 1969 


TO : LEC/teanager, Test Data Reduction Department 672-10 
FROM ED5/Chief , Data Processing Branch 

SUBJECT: Requirements for modification of the THNSKN (Q614) program 


Programing support is required to modify the Thin Skin Transient Heating 
Rate Model Program. Equations 4 and 5 are to be changed to the form shown 
below: 


Equation 4* 


qj = 24kd 



Equation 5?- 



L. 

Mrs, Yvonne Chempinski has started work on these changes in response, to a 
verbal request from Dr. Scott. Work on this task should be a continuing 
charge to project number 4556. For further information, please contact 
E, 0. Grice, extension 5533. ^ 



cc: 

ES56:Dr. C. D, Scott 

ED57/E0Grice (LEC : JMRandal) : be 1 1-20-69 




MUITlN AE MSP ACE SYSTEMS 

16811 EL CAMINO REAL 
HOUSTON, TEXAS 77068 
TELEPHONE i (AREA CODE 713) 4M-00R0 

. LOCKHEED ELECTRONICS COMPANY 

In reply refer to: 
TDR#2272 


December 1, 19 69 


TO: ED5/Chief, Data Processing Branch 

PROM: LEC/Manager, Test Data Reduction Department 

SUBJECT: Technical Memorandum 11-69-636, November 2*1, 1969 

Requirements for Modification of the THNSKN Program 
(Q614) 


Subject Technical Memorandum was received by the Test Data Reduction 
Department on November 28, 1969. 

Programming support for the subject program has been assigned to 
Mrs. Yvonne Chempinski of the Data Acquisition Systems programming 
group. Mrs. Chempinski has completed the requested program modifi- 
cations having worked from a verbal request from Doctor Carl Scott. 

Please direct any additional communications ’ regarding this project 
to Mrs. Chempinski, HU8-0080, extension 345, or to Mr. Len Martin, 
the Engineering Applications Programming team leader, at extension 
336. 

As directed, all work expended on this task will continue to be 
charged to Project Number 4556. 




DL J : JDG : pm 
Distribution: 


\\ 


J. L. Fisher/ED57 
E. 0. Grice/ED57 
J. M. Randal/C 18 



L. Martin/Bll 


J. D. Garrett/Bll 
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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

<>♦** ‘ 

<• V. V* _ \ , 

- • 1 . , 


LYNRON B. JOHNSON SPACE CENTER 

ST ( f u- • * % 

9 %: lj 

i : ' ■ 


Houston, Texas 77053 
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o 

\ v . .■» » 



V** 



^e-ini' 0 

REPLY TO 
ATTN OF: 

r 

U V.* 

December 23, 1974 



Lockheed Electronics Co. Inc . 
■Attn: Mr. Sanies Brose 

Mn.il Stop C-.30 
ISS'il Ei Garni no Real 
Houston, IX 77058 



Dear Mr. Brose: 

A software change request was issued by Mr. E. J. Jung and em on 
October 23, 1974, to incorporate changes into the thin skin heating 
rate data reduction computer program (TIINSKN). Equation (2) of that 
document should be 


24 ' 


rt, N-2 


(T«. - 
10 


r ?: 1 

10 


*T . - T. .) 
U 10' 


The expressions in equations (5) and (6) should be carrying a negative 
sign. 

Please incorporate these corrections into the program. 



C. D. Scott 
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i n r/u h f £ c i i [ c i r c k i c s c a r * :: v. i :: : 

AEROSPACE SYSTEMS DIVISION 

llill (1 (AMINO *1*1 • HOUSTON. MISS Hill • TilUHONC («•(« COOT 111) 411 CM 

In reply refer to: 
#121-8315-199 
J.O. #83-157 

# 

December 31, 1974 

TO: FD5/Chief, Data Processing Branch 

FROM: LEC/Manager, Data Processing Systems Department 

SUBJECT: Transmittal of EXEC VIII Version of the THNSKN Program 


Transmitted herewith is the EXEC VIII version of THNSKN program 
for the Structures and Mechanics System. This version replaces 
the EXEC II version. 


Results from the validation runs, after being checked for 
correctness with the use of LEC generated test data, were sub- 
mitted to the customer, Dr. Carl Scott. Upon his acceptance, 
the Software Change Request was signed by the SIC and NASA 
Project Leader. The Software Change Request form is attached 
giving the descriptive reason for program modification. 


W. fa. Fitzpatrick 



WMF: JEM: jib 
Attachments 


cc: H. Stephenson/FD53 

D. L. McCormick/FD52 
C. Scott/ES3 

E. J. Jung/ES3 

C. E. Hu tchinson/C50 
C. F. Ivea/Cll 


J. Brosc/C30 

J. E. McFarland/C 30 

J. 0. File 


8-23 


I 

i- 


A SUBSIDIARY OF 


0 K H E e 0 A I R C R A F 


C 0 ii P 0 R A T I 0 H 





r-v e-l 


L u C'K H ?. i u 


-> n t f n 

if if i,1 l’ 


h v i 

ft is i . ' 


HIM Cl CAMINO REAL 


ELECTRO t! i C S 

AEROSPACE SYSTEMS DIVISION 

HOUSTON, TEXAS 7 74 SI • TELETHON E (AREA CODE 7U) III 


In reply refer to: 
#121-8315-199 
J.O. #83-157 


December 31, 1974 


TO: FD5/Chief, Data Processing Branch 

FROM: LEC/Manager, Data Processing Systems Department 

SUBJECT: Transmittal of EXEC VIII Version of the THNSKN Program 


Transmitted herewith is the EXEC VIII version of THNSKN program 
for the Structures and Mechanics System. This version replaces 
the EXEC II version. 

Results from the validation runs, after being checked for 
correctness with the use of LEC generated test data, were sub- 
mitted to the customer. Dr. Carl Scott. Upon his acceptance, 
the Software Change Request was signed by the SIC and NASA 
Project Leader. The Software Change Request form is attached 
giving the descriptive reason for program modification. 


I 




W. *N. Fif^patrick 


WMF: JEM: jib 
Attachments 

cc: H. Stephenson/FD53 

D. L. McCormick/FD52 
C. Scott/ES3 

E. J. Jung/ES3 

C. E. Hutchinson/C30 
C. F . Iven/Cll 


J. Brosc/C30 

J. E. McFarland/C30 

J. 0. File 
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Meed Date 


SOFTWARE CHANGE RE CUE ST 


To Bo Completed By Responsible Organizati 
Job Order 83-15 SCR Xn. S3- IS 30 

Program Name THNSKN 

Program No. Q614 

nr* 4300 


Project No mm 

Target Date 3 1/ 7 4 . 


COMPLETE 


D BY RESPONSIBLE ORGANIZATION 


ESTIMATED RESOURCES 


Man-Hours 


Computer Hours Zmi 


REASON TOR CHANGE: 


To improve the accuracy and automation of the computer program THNSKN 
used to calculate heat fluxes from temperature time histories 
measured on thin skin heat transfer models tested in JSC arc tunnel 
facilities . 


CHANGE DESCRIPTION (attach extra sheet if necessary): 
See the attached description. 


iTT.U, TO PRODUCTION: 

ed by:, - 

iTclyrriiATr 

i on approv'd by: /) 

. Li I s-\ / 

C ' * T*-'- ;■ : 

■ n tograt i on Coovdinr tor 


•»\, Project Leader 


LOpnating system 

—New program decks have been sclu titled to production 

Transmittal tape replaces current system CUN Pl’R 

Overlay production CUR/PCK with transmittal CCK/I 

Symbolic elements are contained so file number 

Relocatable elements are contained in file number 

._No changes arc required in execution deck setup 

Changes in run inutntctions are attached 

. -Transmitted prorrain tape number 









A 

. r 


SOFTWARE CHANGE REQUEST 

* l -iJ 


Drrtp Oct, 23, ,1974 


Requester/ • / J / n 

Orgnni 7.nt inn .E.J . Jung/t) . 5 CO t-t / J SC~E S 


Need Date P ec • — . TL- t — I .? 7 . 4 . 


To Be Completed By Responsible Organization 
Job Order 8 3-15 ,SCR Mn - 85-15-50 . . 

Program Name THNSKN 

Program No. Q614 

Project No- 


4300 


Target Date P 2// . 51 /- 7 -l 


REASON FOR CHANGE: 

To improve the accuracy and automation of the computer program THNSKN 
used to calculate heat fluxes from temperature time histories 
measured on thin skin heat transfer models tested in JSC arc tunnel 
facilities. 


CHANGE DESCRIPTION (attach extra sheet if necessary) 


See the attached description. 


TO BE COMPLETED BY RESPONSIBLE ORGANIZATION 



Man-Ilours 1.2 .Q 

Computer Hours 2 . 0 


/■ <sr>i*ech rri Moni toy — 


^Branch Chief 


✓ u o u w 

/s/sZ/v <-/ 


~ 


j ntANSM 1 TTAL TO PRODUCTION : 
Submitted by:, ? 


Pro grammar 


Date 


i '/a-Hda t. i on approved by: /) , 

System n tegrati on Coordinator Date 
. . ‘ »/ # , 
>v / * V • . - / ' • 


_d..Op-5vating system 

New program elects have been sv.bniUed to production 

: Transmittal tape replaces current system CUR, PUR 

Overlay production CUR/PUR wish transmittal CUR/ 1 

Symbolic elements nre contained in file number 

Relocatable elements are contained in file number 

No changes are required in execution deck setup 


!’ 


/ 


.Changes in run instructions are attached 


A, Project Leader 


Date 


Transmitted program tape number 


ACCEPTANCE BY PRODUCTION: 


.* Date to ” t 

ru" Cv.-i -.nleted. 
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Project Leader 


Date 


THIN SKIN PROGRAM MODIFICATION 


Code modifications are required to improve the accuracy and automation 
of the computer program 7IINSKN used to calculate heat fluxes from 
temperature time histories measured on thin skin heat transfer models 
tested in the JSC arc tunnel facilities. The modification described 
here will also allow one discontinuity in the skin thickness (in one 
direction only) and uses finite difference expressions which are 
approximately second order accurate {except at end points) in both the' 
space and time derivatives of the heat conduction equation. 

The following modifications are required: 

1) At present the computer program calculates heat fluxes only 
during one time period. Modify the program so that the number 
of time periods that heat fluxes are calculated during a run 
fascia may be specified in the input. 

2) Co input • the average ‘heat fluxes for each time period and print 
the ?v .rage for each location. 

3) Remove the restriction that if bad thermocouples are at end 
points of rows or columns of thermocouples, the computation is 

or-'-v- o.S 

' aiisorbech Use the next good thermocouple -on the end point. 

4) Compute the conduction correction at end points (last good 

thermocouple) using the same conduction correction as at points 
adjacent to the end point, e.g., ^ ~ ^ 

5) Allowance for discontinuity in skin thickness— In the input 
two skin thicknesses, d] and dz» are t0 l )e specified along with 
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the position y c j of the junction between the two thicknesses 

d = d, for y < y d and d = d 2 for y > y d 

Use the appropriate thickness d for each location. 

6) Compute the mean heat flux for each time interval using the 
formula 



To compute <?i/ at thermocouple locations adjacent to the skin 

'W, ■ 

thickness jump (y d ) use the following formulas for y^Yd^k+l 
ins tea: (4) . 





I d 1 1%. i 

\ -U ->• hu> r 


«►*»*** *»* l ^<* . 

■** | 

w i 
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The bars over the T * s denote time average of the temperatures 
throughout a given interval. The superscripts on the T's in 
the first term denote the first, second, last, and next to 
last temperatures in the time interval. The number of samples 
in a- time interval is N. The first time in the interval is t, 
and the last is tv. The latter two terms in equations (3) and 
(4) are the lateral heat loss terms as used previously except 
now they are averages over the time interval. This procedure 
will replace the one currently used in the code. 


7) Output - Print the start and stop times for each time interval. 
Print out the number of times N in each interval, the times t rt * 


ahH ■hvamrtpv'p H i pqc T at- oarh finw Mn lohnop nv*inf fho boat 

flux terms at each time since these are no longer calculated. 

T • »• X 

Print each of the terms in equation (1): q o r 

bl j -1 S b 


and the net heat flux <?. > . Conti. *i a to compute and print the 

‘J 

linear least squares heat flux in the revised program, Finally 
print the mean of c ^., averaged over the number of time interval 
specified in the input. 

8) Aimiend program documentation to reflect those and any previous 

changes. 


ATTACHMENT 9 


DIAGNOSTIC MESSAGES 


DIAGNOSTIC MESSAGES 



*** FEWER POINTS THAN REQUESTED*** 

This message occurs at the end o£ a time slice when the 
number of .points accumulated for least squares heat flux 
computations is less than the requested (NTME) . This is 
not necessarily abnormal. 

***PARITY ERROR RECORD DELETED*** 


An input record has been skipped because of parity errors. 

END OF FILE ENCOUNTERED ON UNIT 8. PROCESSING OF ACCUMULATED 
DATA WILL BE ATTEMPTED. 

Normally, stop times or stop event levels should be chosen 
so that the end-of-file is never encountered. If at least 
three scans of data have been accumulated, processing will 
be done. 


TOO FEW POINTS AVAILABLE FOR CURVE FIT 

Fewer than three data points have been stored so no least 
squares heat flux computations can be done. 

***ERROR - NX OR NY IS BAD 

Either one or both of NX and NY are too large or too small. 
Currently, the maximum size for either is 50. Either NX or 
NY must have a minimum value of 3. 


X ^ 
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Examples: NX 

NY 


1 

3 

Valid 

3 

1 

‘ Valid 

2 

2 

Invalid 

1 

50 

Valid 

51 

51 

Invalid 


***ERROR - NTME IS LARGER THAN NN 


The requested number of points to be used for least squares 
heat flow computations exceeds array dimensions (NN) . Cur- 
rently NN is 50. 

***ERR0R - TOO MANY CHANNELS 

The number of channels specified for the input tape is too 
large for the array dimensions. Currently, the maximum 
number of channels which can be processed is 260. 
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